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Abstract 

For systems of three identical particles in which short-range forces produce shallow 
two-particle bound states, and in particular for the "pion-less" Effective Field Theory of 
Nuclear Physics, I extend and systematise the power-counting of three-body forces to 
all partial waves and orders, including external currents. With low-energy observables 
independent of the details of short-distance dynamics, the typical strength of a three- 
body force is determined from the superficial degree of divergence of the three-body 
diagrams which contain only two-body forces. This naive dimensional analysis must 
be amended as the asymptotic solution to the leading-order Faddeev equation depends 
for large off-shell momenta crucially on the partial wave and spin-combination of the 
system. It is shown by analytic construction to be weaker than expected in most 
channels with angular momentum smaller than 3. This demotes many three-nucleon 
forces to high orders. Observables like the 4 S3-scattering length are less sensitive to 
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three-nucleon forces than guessed. I also comment on the Efimov effect and limit-cycle 
for non-zero angular momentum. 
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1 Introduction 



Three-body forces parameterise the interactions between three particles on scales much 
smaller than what can be resolved by two-body interactions. Traditionally, they were often 
introduced a posteriori to cure discrepancies between experiment and theory, but such an 
approach is of course untenable when data are scarce or absent, predictive power is required, 
or one- or two-body properties are extracted from three-body data. 

The pivotal promise of an Effective Field Theory (EFT) is that it describes all Physics 
below a certain "breakdown-scale" to a given accuracy with the minimal set of parameters, 
and that it hence predicts in a model-independent way the typical strength with which also 
three-body forces enter in observables. This promise is based on a dimension-less, small 
parameter: the typical momentum of the low-energy process in units of the breakdown- 
scale, namely of the scale on which details of the short-range interactions are resolved. It 
allows one to truncate the momentum-expansion of the forces at a given level of accuracy, 
keeping only and all the terms up to a given order, and thus establishes a "power-counting" 
of all forces. One can then estimate a priori the experimental accuracy necessary to dis- 
entangle particular effects like a three-body force in observables. 

The power-counting is to a high degree determined by naive dimensional analysis [1]. 
As low-energy observables must be insensitive to the details of short-distance Physics, they 
are in particular independent of a cut-off A employed to regulate the theory at short dis- 
tances. Typically, a divergence from loop integrations must therefore be cancelled by at 
least one coefficient C(A) in the EFT Lagrangean, which thus also encodes short-distance 
dynamics. This counter-term enters hence at the same order as the first divergence which 
it must absorb. It ensures that the EFT is cut-off independent at each order, and there- 
fore renormalisable and self-consistent. With the running of the counter-term thus deter- 
mined, its initial condition provides an unknown, free parameter which has to be found 
from experiment. Naive dimensional analysis assumes now that the typical size of the 
counter-term is "natural", i.e. at most of the same magnitude as the size of its running: 
C(A) ~ C(2A) ~ C(2A) - C(A). Thus, it guarantees that the EFT contains at a given 
order the minimal number of free parameters which are necessary to render the theory 
renormalisable, and by that also the minimal number of independent low-energy coefficients 
necessary to describe all low-energy phenomenology to a given level of accuracy. 

When all interactions are treated perturbatively, like in Chiral Perturbation Theory in 
the purely mesonic sector, naive dimensional analysis amounts to little more than count- 
ing the mass-dimension of an interaction [1]. I will demonstrate that such reasoning be- 
comes however too simplistic when some interactions in the EFT must be treated non- 
perturbatively. This is the case in Nuclear Physics, and for some systems of Atomic Physics. 

While the separation of scales between low-energy and high-energy degrees of freedom in 
Nuclear Physics makes it an ideal playground for EFT-methods, finding such a power- 
counting proves a difficult goal for few-nucleon systems. To establish a formalism which is 
self-consistent, agrees with nuclear phenomenology and can firmly be rooted in QCD, one 
has to cope with shallow real and virtual few-nucleon bound-states. The deuteron size of 
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« 4.3 fm for example seems not to be connected even to the soft scales of QCD, e.g. the pion 
mass or decay constant. The effective low-energy degrees of freedom and symmetries of QCD 
dictate a unique Lagrangean; but there are conceptually quite different ordering-schemes 
available for the few-nucleon system which lead to different, experimentally falsifiable pre- 
dictions. The na'ivest versions perturb around the free theory and hence cannot accommo- 
date shallow bound-states at all. They are self-consistent, but not consistent with Nature. 
Weinberg [2] proposed to build few-nucleon systems from a nucleon-nucleon potential which 
consists of contact interactions and pion-exchanges constrained by chiral symmetry. The 
interactions in the potential are ordered following naive dimensional analysis as if the theory 
would be perturbative, and the potential is then iterated to produce the unnaturally large 
length scales in the two-nucleon system by fine-tuning between long- and short-distance 
contributions. Few-nucleon interactions are added using the same prescription. Whether 
this approach correctly and self-consistently reproduces QCD at low energies in the three- 
and more-nucleon sector is an open question. 

However, low-lying few-body bound-states also offer the opportunity for a more radical 
approach: For momenta below the pion-mass, the only forces can be taken to be point- 
like two- and more-nucleon interactions. This Nuclear Effective Field Theory with pions 
integrated out (EFT(^f)) is in the two-nucleon system manifestly self-consistent and proves 
- on quite general grounds - to be the correct version of QCD at extremely low energies, 
once fine-tuning is observed, see Refs. [3-5] for recent reviews. A plethora of pivotal physi- 
cal processes which are both interesting in their own right and important for astrophysical 
applications and fundamental questions, e.g. big-bang nucleo-synthesis and static neutron 
properties, were investigated with high accuracy. One obtains usually quite simple, analytic 
results, and most of the coefficients are determined by simple, well-known long-range ob- 
servables. Recently, a manifestly self-consistent power-counting for the three-nucleon forces 
of EFT(^f) in the 2 Si-wave of A^-scattering was established [6-8]. First high-accuracy cal- 
culations also including external sources are now performed [9]. A remarkable phenomenon 
of this channel is that the first three-body force appears already at leading order to sta- 
bilise the wave-function against collapse [10], leading to a new renormalisation-group phe- 
nomenon, the "limit-cycle" [11-13], manifested also by the Efimov effect [11,14]. This can 
also be shown using a subtraction method [15]. It was also confirmed by an analysis of the 
renormalisation-group flow in the position-space version of the problem [16]. 

EFT(|f) is universal in a dual sense: First, its methods can be applied to a host of 
systems in which short-range forces conspire to produce shallow two-particle bound states: 
One example are identical spin-less bosons, found in bound-states of neutral atoms like 
the 4 He-dimer and -trimer which are bound by van-der-Waals- forces, or loss rates in Bose- 
Einstein condensates near Feshbach resonances, see Ref. [17] for a review. Our results 
are thus readily taken over to such systems. Second, any consistent EFT of nucleons and 
pions must reduce to EFT(^) in the extreme low-energy limit. Therefore, lessons learned in 
the latter shed light on the consistent systematisation of the former. As EFTs are model- 
independent, considerably more sophisticated and computationally involved potential-model 
calculations must agree with the predictions of EFT(f) when they reproduce the two- and 
three-body data which are used as input for EFT (ft) to the same level of accuracy. 
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This article confronts the power-counting of three-body forces in any three-particle system 
with large two-particle scattering lengthes and only contact interactions. It is organised as 
follows: In Section 2.1, the necessary foundations are summarised. After establishing the 
superficial degree of divergence of diagrams which contain only two-body forces in Sect 2.2, 
the far-offshell amplitude of the leading-order Faddeev equation in each partial wave is 
determined analytically in Sect. 2.4 together with a short discussion of the Efimov effect 
in non-integer partial waves. Section 2.4 then classifies at which order a given three-body 
force is needed to render cut-off independent results. Physically relevant consequences are 
discussed in Sect. 3, together with some caveats. After the Conclusions, an Appendix 
sketches some mathematical details. I also correct some errors in a brief summary of some 
of the results in Ref. [18]. 



2 Three-Body Forces in EFT(yf) 

2.1 Three-Body Systems with Large Two-Body Scattering-Length 

We consider three identical particles N of mass M interacting only with contact forces such 
that two particles form a shallow real or virtual two-body bound-state d. As the steps lead- 
ing to the leading-order (LO) scattering-amplitude dN — > dN were often described in the 
literature, they are not covered here; see Ref. [19] also for the notation used in the follow- 
ing. For convenience, a "deuteron" field is introduced as the auxiliary field which describes 
scattering between two particles with an anomalously large scattering length I/7 [20-23]. 
Its propagator is therefore given by the LO-truncation of the effective- range expansion [24]: 

D(q , q) = 1 (2.1) 

7 ~ y t ~ M <7o - ie 

A real bound-state d has at this order the binding energy 7 2 /M <C A^, much smaller than the 
breakdown-scale of EFT on which new degrees of freedom are resolved. In the three-body 
system, an infinite number of diagrams contributes at LO, see Fig. 1. The corresponding 
Faddeev equation for scattering between the auxiliary field d and the remaining particle, 
first derived by Skorniakov and Ter-Martirosian [25], is unitarily equivalent to the original 
problem of scattering between three particles [23,26]. One finds for 7V<i-scattering in the 
/th partial wave in the centre-of-mass system the integral equation for the half off-shell 
amplitude (before wave-function renormalisation) 

00 

tf> (E; k, p) = Stt A K® (E;k,p)--\j 6q q 2 K® (E; q,p) D{E - ^,q) t { ? (E; k, q) . (2.2) 



With the kinematics defined in Fig. 1, E := ffj — jj — ie is the total non-relativistic energy; k 
the relative momentum of the incoming deuteron; p the off-shell momentum of the outgoing 
one, with p = k the on-shell point; and the projected propagator of the exchanged particle 



3 




(ajif *0 (aaf ^ 

Figure 1: Re-summation of the infinite number of LO three-body diagrams (top) into the 
corresponding Faddeev integral equation (bottom). Thick line (D): two-nucleon propagator; 
thin line (/C): propagator of the exchanged nucleon; ellipse: LO half off-shell amplitude. 



on angular momentum I is 

i 



: = \ h n > P '%E + = — Q > (^^) ■ (2-3) 
2 J p z + q z — Mhj + PQ \ PQ J 

-l 

The /th Legendre polynomial of the second kind with complex argument is defined as in [27] 

«w = 5/ * B ■ (2 - 4) 

-1 

The "spin-parameter" A depends on the spins of the three particles and how they combine. 
The values for the physically most relevant three-body systems are summarised in Table 1. 



A = 1 


A = -i 


3 spin-less bosons 


3 nucleons coupled to S = | 


Wigner-symmetric part of 
3 nucleons coupled to S — \ 


Wigner-antisymmetric part of 
3 nucleons coupled to S — » 



Table 1: The spin-parameter A for physical systems of identical particles described by (2.2). 

In EFT(^f), iV<i-scattering in the S = ^-channel is at first sight described by a more com- 
plex integral equation because two-nucleon scattering has two anomalously large scattering 
lengthes: l/^ s in the ^o-channel, and 1/74 in the 3 Si-channel. Therefore, two cluster- 
configurations exist in the three-nucleon system: In one, the spin-triplet auxiliary field d t 
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(the deuteron) combines with the "spectator" nucleoli N to total spin S = | (quartet chan- 
nel) or S = | (doublet channel), depending on whether the deuteron and nucleon spins are 
parallel or anti-parallel. In the other, the spin-singlet auxiliary field d s combines with the 
remaining nucleon to total spin S = |. In the doublet channel, the Faddeev equation is 

thus two-dimensional: The amplitude t® tt stands for the Nd t — > iV<i t -process, rj\ s for the 
Ndt — > iWi s -process, and with D s / t defined analogously to (2.1): 



d,tt 

(0 

d,ts 



j(E;k,p) = 2ir K®(E;k,p) y 



OO 



i J 6qq 2 )C®(E;q,p) 




o 

In the following, we are only interested in the unphysical short-distance behaviour of the 
amplitudes, i.e. in the UV-limit for the half off-shell momenta of (2.5): p,q ^> k, E, 7 a / t . 
This suffices to determine in Sect. 2.4 the order at which divergences need to be can- 
celled by counter-terms parameterising three-nucleon interactions. In this limit, the NN 
scattering-amplitudes are automatically Wignev-SU (4)-symmetric, i.e. symmetric under ar- 
bitrary combined rotations of spin and iso-spin [6, 28, 29]: 

lim D s/t (E- Xj,Q) = lim D(E - Xj,Q) = —7= " ( 2 - 6 ) 

Building the following linear combinations which are symmetric resp. anti-symmetric under 
Wigner-transformations, 

+(Q ._ I ft® _+(Q A +(0 ._ I (*<f) i/(0 A (o 7 \ 

L Wa -~ 2 \ d,tt L d,ts J ' l Wa - — ^ \ d ' U d ' ts ) ' ^ ' 

decouples thus the Faddeev equations of the doublet channel [6]: 
t® \ ( 1 




(p) = 4vr/C«(0;0,p) 

00 

+-^L- |dgg 2 /C«(0;g,p) 


^wa 0De y s m this limit the same integral equation as the quartet-channels, and t® s is identical 
to the one for three spin-less bosons [6, 10]. The problem to construct the UV-behaviour of 
the three-body system with large two-body scattering length simplifies hence to constructing 
the solution of just one integral equation: 

00 

t?(p) = arA taJCW(^- 1 »iM + -j|; (-1)'/ + 1) *?(«) , (2-9) 
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where the amplitude t x (p) depends only on the off-shell momentum p, the partial wave 
Z, and the spin-isospin combination A. The normalisation of the inhomogeneous part only 
provides the overall scale of the solution and is hence irrelevant for the following. 

In a slight abuse of terminology, the names "deuteron" and "nucleon" are used in the 
remainder also when three identical bosons are considered. 

2.2 Divergences and Three-Body Forces at Higher Order 

Naive dimensional analysis is based on the UV-behaviour of the scattering amplitude: As 
outlined in the Introduction, a three-nucleon force is needed at some order as counter-term 
to absorb cut-off dependence induced in the physical amplitudes by divergences, as the 
ingredients of the Faddeev equation are refined to include higher-order effects. The running 
of this three-body force with the cut-off is then assumed to be of the same size as its initial 
condition, which in turn must be determined from a three-body datum. Equivalently, a 
three-body datum is needed at the same order as the first divergence which must be absorbed 
by a three-nucleon interaction. We therefore discuss now the superficial degree of divergence 
of higher-order corrections steming form the "two-body sector" of the theory. 

As will be shown in the next Sub-section, the half off-shell amplitude at large off-shell 
momenta p ^> k is asymptotically given by 

tf\p) oc k l p~ Sl{x) ~ l , (2.10) 

where s/(A) is in general a complex number which depends on the partial wave I and chan- 
nel A. Higher-order corrections to three-nucleon scattering can be obtained by perturbing 
around the LO solution, see Fig. 2. This is numerically tricky [7, 8] also because from 
next-to-next-to-leading order (N 2 LO) on, the full LO-off-shell amplitude must be computed 
and convoluted numerically with the corrections, see the centre bottom graph in Fig. 2. 
However, it allows a simple determination of the order at which the first divergence occurs. 




Figure 2: Top: generic loop correction (rectangle) to the A^d-scattering amplitude at N n LO, 
proportional to q n . Bottom, left to right: Exemplary higher-order contributions to Nd 
scattering from the effective-range expansion (blob) and SD-mixing (circle). Hatched ellipse: 
full off-shell amplitude. Notice that the external legs are on-shell. 
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The asymptotic form of the amplitude at higher orders, and thus its superficial degree of 
divergence, follows from a simple power-counting argument: With q the loop-momentum, 
the non-relativistic nucleon-propagator scales asymptotically as 1/q 2 , its non-relativistic 
kinetic energy as q 2 /M, and a loop integral counts as q 5 /M. The deuteron propagator 
(2.1) approaches 1/q. Only corrections at N n LO which are proportional to positive powers 
of q are relevant in the UV-limit - all other corrections do not modify the leading-order 
asymptotics. They appear together with some coefficients C which encode short-distance 
phenomena and whose magnitude is hence set by the breakdown-scale of the theory. Such 
corrections scale asymptotically as {q/ Kf) n ~ Q n for dimensional reasons. The asymptotics 
of the generic N n LO-correction to the LO amplitude represented by the rectangle in the top 
graph of Fig. 2 is thus proportional to 

k i q -Mx)+i) x q^Ml /q\ n ^ fcJ , ? _ (Si , (y)+1) ^ ^ ^ q n- Sl (x)- Sl ,(X) (2n) 



M q 2 q \A f 

We therefore identify n — S;(A) — sz'(A') as the superficial degree of divergence of a diagram. 
A correction at N n LO diverges when 

Re[n- Sl (X) - s v (\')]>0 . (2.12) 

While S/(A) will turn out to be generically complex, this condition depends only on its 
real part. The power n of the higher-order insertion is on the other hand a positive integer. 
Notice also that this formula is not limited to scattering of three nucleons - it applies equally 
well when external currents couple to nucleons inside the box of Fig. 2. 

Usually, higher-order corrections mix different partial waves, and also Wigner-symmetric and 
Wigner-antisymmetric amplitudes in the doublet-channel. They come from any combination 
of the following effects, some of which are depicted in the lower panel of Fig. 2: 

(1) Effective-range corrections to the deuteron propagator, 



D(qo,q) 



,7- 



7 /<L_ M g 



E 



a (Mq - £) + £ r n (Mq - 

n=l 



-i m 



m=o 



. (2.13) 



With the coefficients r n ~ 1/A^ n+1 of natural size, these contributions are ordered by 
powers of Q ~ q ~ ^Mq : The effective range r enters at NLO as one insertion into 
the scattering-amplitude; Tq at N n LO as n insertions, etc. The correction r n starts 
contributing with one insertion at N 2n+1 LO. They modify the UV-limit of D from 1/q 
in (2.6) at N n LO to q n ~\ 

(2) As two-nucleon forces are non-central, different two-nucleon partial waves mix, e.g. the 
3 Si- and 3 Dx-waves. This leads to mixing and spitting of partial waves also in the 
three-body problem, so that in general S/(A) ^ Sr(A'). By parity-conservation, the 
lowest-order mixing appears for I' — I ± 2. 
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(3) The local vertex of two nucleons scattering via higher partial waves L > contains 
2L positive powers of q and mixes partial waves as well. 



(4) Insertions into D which correct for the explicit breaking of Wignev-SU (4)-symmetry 
are proportional to the differences in the effective- range coefficients of the 1 S - and 
3 Si-system, e.g. to 7 S — jt OT <? 2 (Po,s — Po,t), see also [8]. This leads to mixtures with 
I — I' but A ^ A'. 

(5) Other corrections of less importance, like relativistic corrections to the deuteron prop- 
agator D and to the nucleon propagator. 



A longer remark is appropriate for corrections in which the LO full off-shell amplitude 
is sandwiched between two loop-momenta q±,q2- They lead to overlapping divergences 
in these two variables. Examples are given in the centre bottom graph of Fig. 2, or for a 
correction at N ni+ni LO involving one full off-shell amplitude in Fig. 3. Its asymptotics is 




( 2M ' fc ) 



N ni LO (dr>-«&) (sfr.-ft) N n2 LO 




(l>,\>) 



( 2M ' fc ) 



Figure 3: Exemplary graph containing the full off-shell amplitude. Its kinematics is defined 
after the integrations over the energy-variables gi and g 2 are performed. 



for each off-shell momentum determined by the same Faddeev equation (2.2) with the only 
difference that the cm-energy E and -momentum k become independent variables, with the 
on-shell point at k = p = ^JA{ME + 7 2 )/3. In the asymptotic region E <^ k,p, the integral 
equations for both off-shell momenta in the kinematics defined in Fig. 3 are hence in analogy 
to (2.9) given by: 

oo 




(2.14) 

kp \k pj v^tt J k \q k / 



o 



It is symmetric under the interchange of k and p. In analogy to the solution of the half 
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off-shell amplitude in the next Sub-section, its asymptotics is constructed in App. A. 3 as: 



4\k,p) oc { 



1 




|«l(A) 


kp 


KP, 










kp 


(!) 





for p > k 
for p < fc 



kp 



(2.15) 



The latter follows by analytic continuation to p ~ ^> 7. One may motivate this result 
by observing that in the absence of any other scale, this is the result with both the correct 
mass-dimensions and symmetry-properties. It is independent of the angular momentum and 
spin-parameter. The superficial degree of divergence of Fig. 3 is now easily determined for 
9 ~ qi ~ q 2 scaling alike: 



l.1 ?i M 1 

qi M q\ Ql 



X 



9i 



"1 



x 



(3* 

M ql q 2 



"2 



k l k 1 ' q ni+n2 ~ s '( x )~ s i'( x ') 



x k V q~ MX ' )+1) 



(2.16) 



The overlapping divergence E, k <C 91,92 is hence also included in the previous estimate 
(2.12) when n = n\ + n 2 . One readily generalises to any number j of insertions of full 
off-shell amplitudes with two-nucleon interactions at N ni LO, % = 1, . . . ,j + 1, between them 
and the initial and final half off-shell amplitudes. They all are covered by (2.12), with the 



higher-order correction at N n LO, 



n 



E 



None of the corrections listed above leads at a given order N n LO to stronger modifications of 
the short-distance asymptotics than those induced by effective-range corrections entering at 
the same order, and every order contains also contributions from effective-range corrections. 



2.3 Short-Distance Asymptotics of the Amplitude 

We now just have to determine the unphysical short-distance behaviour of the amplitude to 
infer from (2.12) at which order the first three-body forces are needed to absorb divergences. 
Naively, tf\p) should have the same asymptotics as each of the individual diagrams which 
need to be summed at LO, see top row of Fig. 1. That means, the asymptotic form should 
be given by the inhomogeneous or driving term as 

m 3fc 2 k l 

ty(p) oc lrni/C w (— ; k,p) oc , i.e. s/ iSimp i istic (A) =1 + 1 . (2.17) 

This "simplistic" application of a naive dimensional estimate reflects the expectation that 
three-body forces should enter only at high orders, and that the asymptotics in higher partial 
waves should be suppressed by a centrifugal barrier. Indeed, this estimate would lead from 
(2.12) to the finding that the three-body force in the Zth partial wave - containing at least 21 
derivatives - occurs only at N 2Z+2 LO and is in particular independent of the spin-parameter 
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A. However, the three-body problem consists already at leading order of an infinite number 
of graphs, see Fig. 1. As is well-explored for S-waves, this modifies the solution drastically. 

The integral equation (2.9) can be solved exactly by a Mellin transformation since its 
homogeneous term is scale-invariant and inversion-symmetric; see Appendix A for details. 
An implicit, transcendental, algebraic equation determines the asymptotic exponent s/(A): 



2 i-, A r p±|H] r 



l+s+l I 



1 21 



3 1 
_; 4 



(2.18) 



It depends only on A and I. The function 2 -Pi [a, b; c; x] is the hyper-geometric series [27]. This 
formula comprises the main mathematical result of this article, extends Danilov's result for 
I = [10], and forms in particular the base to power-count all three-body forces. However, 
not all of its solutions solve also the integral equation: While both s and — s are together with 
their complex conjugates solutions to the algebraic equation, only those amplitudes which 
converge for p — > 00 and for which the Mellin transformation exists are permitted. Most 
notably, this constrains Re[s] > —1, Re[s] ^ Re[Z] ± 2; see Appendix A. Furthermore, out 
of the infinitely many, in general complex solutions for given I and A, only the one survives 
as relevant in the UV-limit whose real part is closest to —1, i.e. for which Re[s^(A) + 1] is 
minimal. We consider in the following only those solutions which match these criteria. 

Plots of one of the values in the quadruplet of two-parameter functions {±sj(A), ±s^(A)} 
at fixed I and fixed A, respectively, are given in Figs. 4 and 5. Table 2 lists the first S;(A) 
for the partial waves / < 4 and A = { — |; 1}, compared to the simplistic estimate (2.17). 



partial wave I 


si(A = l) 


*i(A = -§) 


simplistic Z + 1 





1.00624. . . i 


2.16622... 


1 


1 


2.86380... 


1.77272... 


2 


2 


2.82334... 


3.10498... 


3 


3 


4.09040. . . 


3.95931... 


4 


4 


4.96386. . . 


5.01900... 


5 



Table 2: Solutions S;(A) to (2.18) for the most relevant physical systems. 



Let us for the remainder of this sub-section investigate the rich structure of this result. 
Branch points occur e.g. for (/ = 0; A ~ — |), where imaginary parts open. Avoided crossings 
are found e.g. for (I — 1; A = 1), etc. While the solution is in general complex, it is real for 
non-negative integer I in the physical channels discussed above, where A = {1; — =}. The 
only exception is the imaginary solution for (I = 0; A = 1) first found by Danilov [10]. It 
makes a three-body force in this channel mandatory already at LO as the system would 
otherwise be unstable against collapse of its wave-function to the origin, a phenomenon 
well-known to be related to the Thomas- and Efimov-effects [14, 30] and giving rise to a 
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-10 12 3 4 

1 



-10 12 3 4 

1 



Figure 4: The first two solutions S/(A) at A = 1 (left) and A = — |. Solid (dotted): real 
(imaginary) part; dashed: simplistic estimate (2.17); cross (square): real (imaginary) part 
of the asymptotics obtained by a fit of the full solution to the Faddeev equation (2.2) at 
large off-shell momenta. Dark/light: first/second solution. An Efimov effect occurs only for 
|Re[s] | < Re[/ + 1], i.e. when the solid line lies below the dashed one, and Im[s] ^ 0. 
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limit-cycle [11-13] manifesting itself in the Phillips line [6,31,32]. Its interpretation is not 
the scope of this presentation; we only note that the power-counting of three-body forces in 
this channel states that a new, independent three-body force with 21 derivatives enters at 
N 2i LO [8]. It must be seen as coincidence that the naive dimensional estimate in (2.12) - 
where n > was assumed explicitly - leads to the same conclusion. 

In general, s can be complex, as for example at (/ = 0; A < — |), (/ = 1; A > 1) or 
I non-integer, A = {1; — |}. In that case, out of the four independent solutions, the ones 
with Re[s] < —1 must be eliminated as t{p) does not converge for them. In cases like 
(/ G [—0.5819 . . . ; 0.3446 . . .]; A = 1) where all four complex solutions obey Re[s] > —1, only 
the solutions with minimal Re[s;(A) + 1] survive, as shown above. These remaining two 
solutions s := s R ± isj are equally strong and must be super-imposed: 

«?« - si " |s y e-19) 

Usually, Fredholm's alternative forbids that both the homogeneous and inhomogeneous 
integral equations have simultaneous solutions. Therefore, the boundary conditions of the 
integral equation fix the phase S to a unique value. However, when the kernel of the Faddeev 
equation (2.2) is singular, this operator has no inverse and Fredholm's theorem does not 
apply. For this to occur, the solution to the integral equation must be unique only up to a 
zero-mode of the homogeneous version, i.e. 

oo 

8A '■!)' f*rQ>(* + 0'%) = <$<P) (2-20) 
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must have a non-trivial solution a^\q) ^ 0. As the explicit construction in App. A. 2 
demonstrates, this is the case if and only if Re[l + 1] > |Re[s]|, i.e. when Re[s] is smaller in 
magnitude than the blind estimate Si iS j mp ii St ic, (2.17). In that case, a one-parameter family of 
solutions arises. A three-body force is then necessary not to cure divergences but to absorb 
the dependence on the free parameter 5. Its initial condition is not constrained by two- 
body physics but must be determined by a three-body datum. Thus, one finds a limit-cycle 
for such systems, like for three spin-less bosons or the Wigner-symmetric part of the 2 Si- 
wave amplitude in iVci-scattering, (I — 0; A = 1), discussed above. As easily read-up from 
Figs. 4 and 5, this phenomenon occurs for non-negative integer angular momentum only 
when / = and A > 3\/3/(47r), where Re[s] = 0. However, an Efimov effect with complex s 
is often found for non-integer /, e.g. for I G [—0.3544 . . . ; 0.5452 . . . ] in the three-boson case, 
A = 1. A closer investigation will be interesting in view of a conjecture on regularising the 
three-body system in Sect. 3.2. 

A numerical investigation of the Faddeev equation (2.2) confirms these findings. In order 
to compute a solution, one introduces a cut-off A which is un-physical and thus not to be 
confused with the breakdown-scale of the EFT. The numerical values of sz(A) are found 
from fitting the half off-shell amplitude t®(E, k;p) at E, ^, 7 -C p C A to the asymptotic 
forms (2.10) and (2.19). A grid of 100 points is easily enough for a numerical precision in s 
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of about 1% 1 . Agreement between the numerical and analytical solution is excellent also 
at non-integer I and A ^ {1, — |}, see Fig. 6 besides Figs. 4 and 5 for examples. Particularly 
interesting is in that context the neighbourhood around the A = — ^-solution in the S-wave 
channel, / = 0. Here, the first solution to the algebraic equation (2.18) is s = 2, but the 
Mellin transform does not exist at that point because Re[s] = I ±2. The system is here also 
close to the branch-point at (A = —0.50416 . . . ; s ~ 2.0836 . . . ), where an imaginary part 
opens for smaller A. Another branch-point lies at (/ = 1; A = 1.0053 . . . ; s = 2.93164 . . . ). 




Figure 6: Left: Numerical and analytical solution for sj(A) at I — around A = — \. Right: 
Numerical determination of S;(A), exemplified for / = 1, A = 1.5, comparing "data" (crosses) 
and the fitted function (2.19) (solid line). Notation as in Fig. 4. 



To summarise, the algebraic equation (2.18) for s;(A) provides asymptotic solutions of the 
form (2.10) to the three-body Faddeev equation (2.2) for Re[s] > —1 and Re[s] ^ Re[/±2]. 
Only those solutions are relevant in the UV-limit p S> 7, E, k for which Re[s + 1] is minimal. 
The Efimov effect occurs only if Im[s] 7^ and | Re [s] | < Re[/ + 1], because only then is the 
kernel of the integral equation not compact. 

2.4 Ordering Three-Body Forces 

Although divergences can occur as soon as the superficial degree of divergence Re[n — 
s/(A) — si'(X')] > 0, only those are physically meaningful which can be absorbed by three- 
body counter-terms, i.e. by a local interaction between three nucleons in the given channels. 
Naive dimensional analysis does not construct the three-body forces. It thus predicts some 
divergences which are absent when the diagram is actually calculated. There is for example 
no Wigner-S , f/(4) anti-symmetric three-body force without derivatives [6], so that the diver- 
gence must in this case be at least quadratic for infinite cut-off, Re[n — si(X) — sj'(A')] > 2. 
For finite cut-off A, a three-body force without derivatives can be constructed which is 
non-local on a scale smaller than 1/A - but appears of course local at scales smaller than 

1 A simple Mathematica-code can be down-loaded from http://www.physik.tu-muenchen.de/~hgrie. 
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the break-down scale < A of EFT(^f). As its coefficient dis- appears when the cut-off is 
sent to infinity, it is in a renormalisation-group analysis classified as "irrelevant". Except for 
this example which is relevant below, this article will however simply assume that the first 
three-body force enters with the first divergence. To construct these forces in detail is left 
for future, more thorough investigations. 

The solution of (2.18) approaches for large integer / in the physically most interesting cases 
A = {1; — |} the simplistic estimate of (2.17): s s i mp n st i c = 1 + 1; see Table 1 and Figs. 4, 
5. This reflects that the Faddeev equation should be saturated by the Born approxima- 
tion in the higher partial waves because of the ever-stronger centrifugal barrier between the 
deuteron and the nucleon. Therefore, three-body forces enter in most channels for all prac- 
tical purposes at the same order as suggested by the simplistic argument, namely +2 LO 
between the Ith and I'th partial waves. It is therefore convenient to introduce the variable 

A,(A) :=s,(A)-(Z+l) (2.21) 

which parameterises how strongly simplistic and actual asymptotic form differ. For A > 0, 
the superficial degree of divergence of the LO amplitude is weaker than guessed by (2.17). 

In the lower partial waves I, I' < 2, however, the blind expectation deviates substantially 
from the exact solution; see Table 3. For (/ = 0; A = 1), for example, sz, S i mp iistic = 1+1 under- 
estimates the short-distance asymptotics of t(p), while s (l) = 1.006 ... i is even imaginary. 
A limit-cycle signals that one must include a three-body force already at LO, as briefly 
hinted upon above. Multiple insertions of three-body forces are not suppressed. 

For two partial waves, the formula (2.17) substantially over-estimates the asymptotics 
oftf\ Therefore, a three-body force is in channels which involve these partial waves weaker 
than predicted by a simplistic application of naive dimensional analysis. Consider first the 
case of three bosons with (I = 1; A = 1): s = 2.86 . . . > s S j mp ij St ic = 2, A = 0.86 .... While 
the first divergence from the two-body sector arises in this partial wave at N 5 72 LO, one 
would - following (2.17) - have predicted the first three-body force as necessary already at 
N 4 LO. It is in this channel hence demoted by ~ 1.7 orders. 

The situation is even more drastic in the 4 S3-channel of A^-scattering, (/ = 0; A = — |): 
Here, only divergences which are at least quadratic are physical because the first three-body 
force must contain at least two derivatives since the Pauli principle forbids a momentum- 
independent three-nucleon force - or equivalently, no Wigner-SU (4) anti-symmetric three- 
body force exists [6]. Therefore, the divergence condition (2.12) reads Re[n — si(X) — si'(X') > 
2. Since s = 2.16 . . . > s S im P iistic = 1, the first three-body force enters thus not at N 4 LO but 
at least two orders higher, namely at N 633 ' LO. 

As as example for mixing between partial waves, consider the 4 S3-wave: It mixes with 

both the 4 D|-wave (A = — |) and the Wigner-symmetric and anti-symmetric components 
of the 2 D|-wave (A = 1 or — |). All of them are already close to the estimate s/ iS j mp ii st i c = 
/ + 1 = 3 > s/=o(A = — 5). Still, the first divergences induced by this mixing start from 
(2.12) at N~ 5 LO, i.e. approximately one order higher than blindly guessed. 

More modifications induced by mixing and splitting of partial waves as well as explicit 
breaking of the Wigner-S77(4)-symmetry are straight-forwardly explored, but left to a future 
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channel 


estimate 
Re[s ; (A) + s v (\')] 


simplistic 

I + r + 2 




(A;/) (A';/') 


partial waves 


(1;0) (1;0) 

(i;0) (-|;0) 

(-|;0) (-|;0) 


2q 2q 
2c 2q 
2q 2q 


LO 

N 2 - 2+2 LO 
N 4 - 3+2 LO 


N 2 LO 
N 2+2 LO 


promoted 
demoted 


(i;0) (-|;2) 
(-|;o) (-|;2) 


2 Sws- 4 D 

2 Swa" 4 D 


N 31 LO 
JN°°LO 


N 4 LO 


promoted 
demoted 


(i;i) (-i; 1 ) 
H;i) 


2p 2p 
r Ws- r Ws 

2p 2p 2p 4p 
r Ws- r Wa , r Ws- r 

2p 2p 4p 4p 
r Wa" r Wa i -T- r 


N 5 - 7 LO 
N 4 - 6 LO 
JN LO 


N 4 LO 


demoted 
demoted 


(-|;0) (-|;0) 
(-|;0) (1;2) 

H;0) (-|;2) 


4g_4 S 

4q 2p. 

4 S- 2 D Wa , 4 S- 4 D 


N 43+2 LO 
N 5 -°LO 
N 53 LO 


N 2+2 LO 
N 4 LO 


demoted 
demoted 
demoted 


(i;2) (i;2) 
(i;2) (-|;2) 

(-|;2) (-|;2) 


2 Dws- 2 Dws 

2p. 2p. 2p. 4p. 

2 Dwa- 2 D\va j 4 D- 4 D 


N 5 - 6 LO 
N 5 - 9 LO 
N 62 LO 


N 6 LO 





Table 3: Order n at which the leading three-nucleon force enters for the lowest channels 
I, I' < 2, comparing the simplistic estimate (2.17) and the actual values (2.12/2.18). The list 
follows the physical partial wave mixing, and the sub-script Ws (Wa) denotes the Wigner- 
symmetric (anti-symmetric) contribution. In the 2 Sw a - and 4 S3-channels, the absence of 
a three-nucleon force without derivatives is taken into account by the factor "+2". The 
last column indicates whether the three-body force is stronger ("promoted") or weaker ("de- 
moted") than the simplistic estimate suggests. When the difference between the two is in 
magnitude smaller than 0.5, they are quite arbitrarily assumed to enter at the same order. 

publication. Table 3 summarises the findings for the physically most relevant three-body 
channels I, I' < 2. Except in the 4 S3-channel, it does not take into account whether a three- 
body counter-term can actually be constructed at the order at which the first divergence 
occurs. However, while this can make a three-body force occur at a higher absolute order 
than listed, the relative demotion of a three-body force to higher orders by modifications of 
the superficial degree of divergence holds. 

"Fractional orders" are a generic feature of (2.18), combined with (2.12). Consider again 
as example the case (I = 0; A = — =), where the first two-body divergence appears formally 
at n — 6.33 . . . , while including a two-body correction with fractional order is of course 
impossible: Clearly, the N 7 LO-amplitude diverges without three-body forces, but one could 
also argue that it is prudent to include a three-body force already at N 6 LO because the 
higher-order correction to the amplitude converges only weakly, namely as g -0,33 . Therefore, 
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the integral over q becomes unusually sensitive to the amplitude at large off-shell momenta, 
above the breakdown-scale of the EFT, and therefore to details of Physics at distances 
on which the EFT is not any more valid. In the more drastic case S/ =4 (A = — |) = 5.02 . . . , 
A = 0.02. . . , the first three-body force enters following (2.12) at N 1 LO. Therefore, no 
divergence arises in the two-body sector before N 11 LO, but it seems reasonable to include it 
already at N 10 LO because the higher-order corrections from two-nucleon insertions converge 
then generically to a cut-off independent result only very slowly, namely as g~°' 04 . 

Naive dimensional analysis cannot decide the question at which order precisely a "frac- 
tional divergence" gives rise to a three-body force as it argues on a diagram-by-diagram 
basis, missing possible cancellations between different contributions at the same order. One 
way to settle it is to see whether the cut-off dependence of observables follows the pattern 
required in EFT. Recall that N n LO corrections contribute to observables typically as 

compared to the LO result and that low-energy observables must be independent of an 
arbitrary regulator A up to the order of the expansion. In other words, the physical scat- 
tering amplitude must be dominated by integrations over off-shell momenta q in the region 
in which the EFT is applicable, q < Kf. As argued e.g. by Lepage [33], one can therefore 
estimate sensitivity to short-distance Physics, and hence provide a reasonable error-analysis, 
by employing a momentum cut-off A in the solution of the Faddeev equation and varying 
it between the breakdown-scale A^ and oo. If observables change over this range by "con- 
siderably" more than Q n+1 , a counter-term of order Q n should be added. This method is 
frequently used to check the power-counting and systematic errors in EFT (ft) with three 
nucleons, see e.g. most recently [19]. A similar argument was also developed in the context 
of the EFT "with pions" of Nuclear Physics [34,35]. Such reasoning goes however beyond 
the clear prescription according to which only divergences make the inclusion of counter- 
terms mandatory and opens the way to a softer criterion - which is obviously formulated 
rigorously only with great difficulty. How to treat "fractional orders" in a well-prescribed 
and consistent way must thus be investigated further. 

2.5 How Three-Body Forces Run 

Before turning to practical consequences of these observations, let us for a moment inves- 
tigate how the strengthes of three-body forces have to scale with q in order to absorb the 
divergences (2.11) from two-body interactions. In contradistinction to the above consid- 
erations where the specific form of the three-body force did not enter, we now limit the 
discussion to those three-body forces which can be re-written as deuteron-nucleon interac- 
tions. Clearly, all three-body forces which are needed to absorb divergences from two-nucleon 
effective-range corrections proportional to r n fall into that class 2 . As discussed in Sect. 2.2, 
this is no severe restriction because every divergence contains such a piece. The leading con- 
tributions are given by the diagrams of Fig. 7. At even higher orders, two- and three-body 

2 But not necessarily three-body forces which contribute to the mixing and splitting of partial waves. 
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Figure 7: Generic leading corrections from three-body forces which can be re-written as 
iVd-interactions with at least I + I' derivatives. 



corrections occur simultaneously in one graph. Let the three-body force between deuteron 
and nucleon with relative incoming momentum p and outgoing momentum p' scale as 



PV' ( 2 - 2 3) 

i 

where the dimension-less coupling h encodes the short-distance details of three interacting 
nucleons which are not resolved in EFT (ft). It absorbs hence also the divergences generated 
at order S/(A) + si> from two-body interactions, (2.11), to render the result at this order 
insensitive of unphysical short-distance effects. The parameter h must thus formally scale 
as q to some power a which is determined such that at least one of the three-body force 
graphs appears at the same order as the divergence. The graphs scale and diverge as: 

(a) ~ k l k 1 ' q a , no divergence 

(bl) ~ k l k v q a ~ A ^ , diverges for RefAj(A)] < 

(2.24) 

(b2) ~ k l k 1 ' g"" A i'( A ') , diverges for Re[A,,(A')] < 

(c) ~ k l k 1 ' g «-A i (A)-A ; /(A') ^ diverges for Re[A,(A)] < or Re[A,,(A')] < 

The tree-level contribution is of course free of divergences. Notice that the three-body 
force h absorbs for non-integer s also the non-analytic piece of the divergence (2.11), and 
in particular the phase 5 when Im[s] ^ 0, see (2.19). This piece is non-analytic in the 
un-physical off-shell momentum q, but of course analytic in the low-energy momentum k. 

The graphs containing three-body forces enter at the same order for all channels which 
follow the simplistic estimate of (2.17), i.e. A/(A), A//(A') = 0. Then, all three-body cor- 
rections (a-c) occur at the same order a, the three-body force counts as h ~ q a = q l+l+2 i 
and the logarithmic divergences of the loop-diagrams (bl/2,c) are absorbed into h as well. 
To determine the order at which a three-body force enters, it is therefore sufficient to count 
in this case its mass-dimension, which is also given by I + I' + 2, see (2.23). This is nearly 
realised in the higher partial waves, where |A| is usually not bigger than 0.3. 

For (/ = 0; A = 1), however, Re [A] = —1 and a = 0. Now, multiple insertions of three- 
body forces are not suppressed and the Efimov effect mandates including the three-body 
force in the LO Faddeev equation. The strength h depends on the arbitrary phase 5, showing 
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a limit-cycle [11-13] as manifested in the Phillips line [6, 31, 32]. Interestingly, the diagrams 
(a,bl/2) - as well as their analogues with higher-order three-nucleon forces - now become 
following (2.24) formally corrections of higher order, as found numerically in Ref. [7]. 

As seen in the previous Section, three-body forces appear in many channels at higher 
orders than expected. This now also re-groups the graphs containing three-body forces. 
According to the scaling properties (2.24), the tree-level diagram (a) is in the (Z = 1; A = 1)- 
channel « 1.7 orders weaker than the leading three-body diagram (c) because A = 0.86 
It can therefore safely be neglected when absorbing the leading divergences from two-body 
insertions. The graphs (bl/2) are down by ~ 0.9... orders. For (I — 0; A = — |), this 
is even more pronounced: With A = 1.16..., the tree-level three-body contribution (a) 
is suppressed by more than two, and (bl/2) by more than one order against the sand- 
wiched three-body graph (c), so that both can be neglected when the leading divergences 
are absorbed into (c) only. Notice that all three-body corrections converge for A > 0. 

Possible overlapping divergences complicate a similar analysis in the case of three-body 
forces which cannot be re-written as iV<i-interactions, warranting further investigations 
which are however not central to this presentation. 

3 Consequences 

The first goal of this publication has been reached: Eq. (2.12) is an explicit formula for the 
order at which the first three-body force must be added to absorb divergences. It depends on 
the partial wave / and channel A via the exponent S/(A) which characterises the asymptotic 
form of the half off-shell scattering matrix tf\E,k;p) and is determined by (2.18). A 
simplistic application of naive dimensional analysis, (2.17), provides a good estimate for all 
partial waves I > 2. However, it over-rates three-body forces of the three-nucleon system for 
example in the 4 S3- and 2 P-channels, while it under-estimates them e.g. in the 2 Si-wave; 
see the summary in Table 3. In the case of three spin-less bosons, the P-wave three-body 
interaction is weaker, while the S-wave interaction is stronger than the simplistic argument 
suggests. With these findings, the EFT of three spin-less bosons and the pion-less version 
of EFT in the three-nucleon system, EFT(^f), are self-consistent field theories which contain 
the least number of counter-terms at each order to ensure renormalisability. Each three- 
body counter-term gives rise to one subtraction-constant which must be determined by a 
three-body datum. Let us explore in this Section some physically relevant results which can 
be derived from these findings. 

3.1 Context 

Amending Naive Dimensional Analysis: As outlined in the Introduction, power- 
counting by naive dimensional analysis amounts for perturbative theories to little more 
than counting the mass-dimensions of the interactions [1]. In this case, only a finite number 
of diagrams contributes at each order. When the LO amplitude is however non-perturbative, 
i.e. an infinite number of diagrams must be summed to produce shallow bound-states, then 
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the situation changes: The LO amplitude can follow for large off-shell momenta a different 
power-law than the one which one obtains when one considers the asymptotic form of each 
of the diagrams separately. Then, the "canonical" application turns out to be too simplistic 
and must be modified as in Sect. 2.3. This is neither a failure of naive dimensional analysis, 
nor should it come completely unexpectedly. An example is indeed already found in the 
two-body sector of EFT(^). Recall that an infinite number of two-body scattering diagrams 
are re-summed into the LO deuteron propagator (2.1) to produce the shallow two-body 
bound state, Fig. 8. Each of the diagrams diverges as g n , with n the number of loops. Their 




Figure 8: Re-summation of the infinite number of LO two-body diagrams into the deuteron 
propagator given by (2.1). 



sum converges however as 1/q for large momenta, see (2.6). In this case, the solution is 
obtained by a geometric series and the necessary changes are easily taken into account, see 
the reviews [3-5, 17]. What may come as a surprise is that this can also happen when both 
all LO diagrams separately and their re-summation are ultra-violet finite. In the three-body 
case, all diagrams actually show the same power-law behaviour l/q l+1 , see (2.17). However, 
the re-summed form looks very different, exhibiting a non-integer and even complex power- 
dependence (2.18). This does not occur for every system with shallow bound-states. For 
example, the exchange of a Coulomb photon between two non-relativistic, charged particles 
in non-relativistic QED scales asymptotically as 1/q 2 . The exact solution of the Coulomb 
problem has the same scaling behaviour. 

External currents: The power-counting of three-body forces developed above applies 
equally when external currents couple to the three-nucleon system. The only change is 
that the higher-order interaction in Fig. 2 becomes more involved, introducing also the 
momentum- or energy-transfer from the external source as additional low-energy scales. 

Three-body forces at higher orders: Another trivial extension is to power-count three- 
body forces beyond the leading ones. In that case, the superficial degree of divergence must 
by analyticity be larger than a positive even integer, Re[n — sz(A) — sz'(A')] > 2m, m G 2No. 
The higher-order three-body force contains 2m derivatives more than the leading one (m = 
0) and enters 2m orders higher. We used this already to power-count the first three-body 
force in the 4 S|-channel. The power-counting based on naive dimensional analysis agrees for 
the 2 Si -channel with the one which was recently established by a more careful and explicit 
construction [6-8]. 
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3.2 Conjectures 

Predicting the 4 S|-scattering length: The first conjecture follows from the observation 
above that three-body forces are demoted in the 4 S|-wave from a N 4 LO-effect by two orders 
to N 63 LO. This has immediate consequences for the quartet-S wave scattering-length of the 
nucleon-deuteron system which has drawn substantial interest recently. Its knowledge sets at 
present the experimental uncertainty in an indirect determination of the doublet scattering 
length [36], which in turn is well-known to be sensitive to three-body forces [31]. It was 
determined repeatedly in EFT(^) at N 2 LO, with different methods to compute higher-order 
corrections agreeing within the predicted accuracy [21-23,26], e.g. most recently [19]: 

a( 4 S|) = (5.091 + 1.319 - 0.056) fm = [6.35 ± 0.02] fm . (3.1) 

The theoretical accuracy by neglecting higher-order terms is here estimated conservatively 
^ Q ~ 7 A^~m GV ~ I of the difference between the NLO- and N 2 LO-result. This agrees very 
well with experiment [37], [6.35 ± 0.02] fm, albeit partial wave mixing, iso-spin breaking 
and electro-magnetic effects are not present in EFT(^) at N 2 LO. As the amplitude decays 
at large off-shell momenta as 1/p 316 ", see Table 1, it is not surprising that a( 4 Si) is to a 
very high degree sensitive only to the correct asymptotic tail of the deuteron wave function. 
The first three-body force 3 enters not earlier than N 6 LO - taking a conservative approach 
as described above to round the "fractional order". Indeed, if the theoretical uncertainty 
decreases steadily from order to order as it does from NLO to N 2 LO, then one should be 
able to reach an accuracy of ±(|) 4 x 0.06 fm < ±0.001 fm with only two-nucleon scattering 
data as input - provided those in turn are known with sufficient accuracy. Indeed, this is not 
much smaller than the range over which modern high-precision potential-model calculations 
differ: [6.344. . .6.347] fm [38,39]. To use this number hence as input into a determination 
of the doublet scattering-length as a( 2 Si) = [0.645 ± 0.003(exp) ± 0.007(theor)] fm [36] 
seems justified, and the error induced by the theoretical uncertainty might actually be over- 
estimated. Notice that if the three-nucleon force would occur in EFT(f) at N 4 LO as the 
simplistic expectation (2.17) suggests, the error should be of the order of (|) 2 x 0.06 fm « 
0.007 fm, considerably larger than the spread in the potential-model predictions. 

An alternative regularisation: More speculative is the possibility for a new regularisa- 
tion scheme. In principle, (2.18) gives the asymptotics S/(A) of the half off-shell amplitude 
for arbitrary - even complex - / and A. As this function is largely analytic, one could use 
analytic continuation for a "partial wave regularisation" of the three-body system. This is 
particularly attractive to regulate the limit-cycle problem of the 2 Si-wave (/ = 0; A = 1), 
whose practical implications are at present mostly discussed by cut-off regularisation. How- 
ever, the algebraic solution to (2.18) suffers - as discussed in Sect. 2.3 - from constraints 
by branch-cuts and regions where the Mellin transformation does not exist. In addition, a 

3 In the mixing between the 2 Si-, 2 D- and 4 D-waves, three-body forces appear already at N 5 LO, see 
Table 3, but they are irrelevant for the scattering length. 



20 



limit-cycle is encountered only when s has an imaginary part, and its real part is smaller 
than the simplistic estimate (2.17). However, the imaginary part dis-appears in the vicin- 
ity of the physical point (/ = 0; A = 1) only where s has a branch-cut, see Figs. 4 and 5. 
Blankleider and Gegelia [40] attempted to use analyticity in A at fixed / = to regulate the 
three-boson problem at LO without resorting to three-body forces to stabilise the system. 

3.3 Caveats 

Weak points in the derivation should also be summarised: 

(1) As always in naive dimensional analysis, one obtains only the superficial degree of 
divergence. Except in the 4 S3-channel, the classification of Table 3 does not take into 
account whether a three-body counter-term can actually be constructed at the order 
at which the first divergence occurs. However, while the actual degree of divergence 
must usually be determined by an explicit calculation, it is is never larger than the 
superficial one. Therefore, a three-body force can possibly occur at a higher absolute 
order than predicted by naive dimensional analysis, but this applies then equally well 
to the simplistic estimate. Therefore, the relative demotion of a three-body force to 
higher orders by the modified superficial degree of divergence holds. In this context, 
the three-body forces which can contribute in a given channel, and in particular to 
partial-wave mixing, should be constructed explicitly. Here, the symmetry principles 
invoked above can be helpful. 

(2) The divergence of each diagram was considered separately, missing possible cancel- 
lations between different contributions at a given order. This would again demote 
three-body forces to higher orders than determined by the superficial degree of diver- 
gence. It would also require a fine-tuning whose origin would have to be understood. 
When further re-summations of infinitely many diagrams should be necessary beyond 
LO, naive dimensional analysis must be amended further. This could happen when 
the power-counting developed here does not accord to nuclear phenomenology. 

(3) Modifications by overlapping divergences should also be explored. Again, they weaken 
the degree of divergence, but are particularly important for those three-body forces 
which cannot be re-written as iV<i-interactions. 

(4) The problem of "fractional orders": Since the LO amplitude involves an infinite num- 
ber of graphs, equivalent to the solution of a Faddeev equation, the amplitude ap- 
proaches for large half off-shell momenta generically a power-law behaviour g^'M" 1 
with irrational and even complex powers. Three-body forces contain therefore follow- 
ing Sect. 2.5 non-analytic pieces. We assume that this behaviour is changed at higher 
orders only by integer powers because higher-order corrections involve only a finite 
number of diagrams after the LO-graphs are summed into the Faddeev equation. At 
which concrete order a given three-nucleon interaction needs to be included to ren- 
der observables cut-off independent can therefore become a question beyond the clear 
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prescription according to which only divergences make the inclusion of counter-terms 
mandatory. A softer criterion is formulated rigorously only with great difficulty. 

(5) The two-nucleon propagator (2.1) at the starting point of the derivation was taken to 
be already renormalised. This should pose no problem as the Faddeev equation was 
solved without further cut-offs, so that no overlapping divergences occur. Indeed, any 
"scale-less" regulator (like dimensional regularisation) will lead to the same result. 

(6) Partial wave-mixing and -splitting as well as mixing between Wigner-S77 (4) symmetric 
and anti-symmetric amplitudes should be considered in more detail. However, we saw 
already examples where these effects are suppressed. For example, the power-counting 
in the partial wave with the lowest angular-momentum amongst those that mix is 
unchanged, see the 4 S|- 4 D|- 2 D|-mixing in Table 3. The higher the partial wave, the 
closer is its asymptotics to the simplistic expectation (2.17). 

(7) We assumed - as usual in EFT - that the typical size of three-nucleon counter-terms 
is set by the size of their running. There are cases where the finite part of a counter- 
term is anomalously large and thus should be included already at lower orders than 
the naive dimensional estimate suggests. One example is the anomalous iso-vector 
magnetic moment of the nucleon which is as large as the inverse expansion parameter 
of EFT(^f), Ki = 2.35 ~ 1/Q [41]. Such cases are however rare and must be justified 
with care. 



Finally, Blankleider and Gegelia [40] claimed in an unpublished preprint 5 years ago 
that the 2 Si-wave problem can be solved at LO without resorting to a three-body force to 
stabilise the system against collaps. According to them, if the Faddeev equation has multiple 
solutions, then only one is equivalent to the series of diagrams drawn in Fig. 1. We focus in 
this article on the higher partial waves where the Faddeev equation has - as demonstrated 
in Sect. 2.3 - always unique solutions for integer I > and A G {1; — |}, so that the alleged 
discrepancy cannot arise. Distracting the reader for a moment, one may however point out 
a few observations which contradict the claim of Ref. [40]. The derivation of the Faddeev 
equation is just a special case of Schwinger-Dyson equations, which are well-known to be 
derived in the path-integral formalism without resort to perturbative methods, see e.g. [42, 
Chap. 10]. One resorts to a "series of diagrams" only for illustrative purposes like in Fig. 1. 
Recall that in the case of the three-body system, no small expansion parameter exists in 
which this series can be made to converge absolutely. In addition, and on a less formal level, 
well-known properties of the three-body system like the Thomas and Efimov effects [14, 30] 
and the Phillips line [6,31,32] are not explained under the assertions of Ref. [40]. These 
universal properties were recently also tested experimentally, e.g. for particle-loss rates in 
Bose-Einstein condensates near Feshbach resonances, see [17] for a review. 
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4 Conclusions and Outlook 



In this article, the ordering of three-body contributions in three-body systems also coupled 
to external currents was constructed systematically for any EFT with only contact interac- 
tions and an anomalously large two-body scattering length. Evading explicit calculations, 
the result is based on naive dimensional analysis [1], improved by the observation that be- 
cause the problem is non-perturbative already at leading order, the solution to the Faddeev 
integral equation does for large off-shell momenta not follow a simplistic dimensional es- 
timate, Sects. 2.3 and 3.1. This was shown by constructing the analytical solution to the 
Faddeev equation in that limit for arbitrary angular momentum and spin-parameter. One 
could thus develop a "partial-wave regularisation" as an alternative to regulate and renor- 
malise the three-body system. A simplistic approach to naive dimensional analysis fails for 
systems which are non-perturbative already at leading order. 

In order to keep observables insensitive to the details of short-distance Physics, one 
employs the canonical EFT tenet that a three-body force must be included if and only if it 
is needed as counter-term to cancel divergences which can not be absorbed by renormalising 
two-nucleon interactions. After determining the superficial degree of divergence of a diagram 
which contains only two-nucleon interactions in Sect. 2.2, this was used in Sect. 2.4 to classify 
the relative importance of three-body interactions for each channel, also for partial-wave 
mixing and splitting. With these results, the EFT of three spin-less bosons and EFT(f) 
become self-consistent field theories which contain the minimal number of parameters at each 
order to ensure renormalisability and a manifest power-counting of all forces. Each such 
three-body counter-term gives rise to one subtraction-constant which must be determined 
by a three-body datum. 

It must again be stressed that three-body forces are in EFT(^) added not out of phe- 
nomenological needs. Rather, they cure the arbitrariness in the short-distance behaviour 
of the two-body interactions which would otherwise contaminate the on-shell amplitude, 
and hence make low-energy observables cut-off independent on the level of accuracy of the 
EFT-calculation. Recall that the theory becomes invalid at short distances as processes 
beyond the range of validity of EFT(^f) are resolved, namely the pion-dynamics and quark- 
gluon sub-structure of QCD. Three-body forces are thus not introduced to meet data but 
to guarantee that observables are insensitive to off-shell effects. Only the combination of 
two-body off-shell and three-body effects is physically meaningful. 

Most of the three-nucleon forces in partial waves with angular momentum less than 3 
have a weaker strength than one would expect from a blind application of naive dimen- 
sional analysis, see Table 3. This might seem an academic dis-advantage - to include some 
higher-order corrections which are not accompanied by new divergences does not improve 
the accuracy of the calculation; one only appears to have worked harder than necessary. 
However, it becomes a pivotal point when one hunts after three-body forces in observables: 
In order to predict the experimental precision necessary to dis-entangle these effects, the 
error-estimate of EFT is a crucial tool. For many problems, this makes soon a major differ- 
ence in the question whether an experiment to determine three-body force effects is feasible 
at all. One such consequence was discussed in Sect. 3.2: The 4 S3-wave scattering length is 
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fully determined by two-nucleon scattering observables on the level of ±0.001 fm according 
to the power-counting of EFT(f) developed here. This is more than a factor of ten more 
accurate than the present experimental number [37], but supported by the observation that 
all modern high-precision two-nucleon potentials predict this observable to a similarly high 
accuracy [38, 39]. If a three-nucleon force (even one saturated by pion-exchanges) would oc- 
cur at the order at which it is blindly expected, the spread in the potential-model predictions 
should be considerably larger. 

The consequences for other observables like the famed Ay-problem [43] should also be 
explored. This will be particularly simple in EFT(f) because the theory is less involved and 
Table 3 sorts the three-nucleon forces according to their strengthes, indicating also their 
symmetries and the channels in which they contribute on the necessary level of accuracy. 
The conclusions, conjectures and caveats of Sect. 3 summarise a number of further interesting 
directions for future research. 

To classify the order at which a given two- or three-nucleon interaction should be added in 
Chiral EFT, the EFT of Nuclear Physics with pions as explicit degrees of freedom, I suggest 
to follow a path as for EFT(^) which complements the so-far mostly pursued phenomeno- 
logical approach: At leading order, the theory must be non-perturbative to accommodate 
the finely-tuned real and virtual two-body bound states in the S-waves of two-nucleon scat- 
tering. After that, only those local two- and three-nucleon forces are added at each order 
which are necessary as counter-terms to cancel divergences of the amplitudes at short dis- 
tances. This mandates a more careful look at the leading-order, non-perturbative scattering 
amplitudes to determine their ultraviolet-behaviour and superficial degree of divergence, 
see e.g. [44] and references therein. It leads at each order and to the prescribed level of 
accuracy to a cut-off independent theory with the smallest number of experimental input- 
parameters. The power-counting is thus not constructed by educated guesswork but by 
rigorous investigations of the renormalisation-group properties of couplings and observables 
by EFT-methods. Work in this direction is under way, see also [45], and the future will 
show its viability. 
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A Solving the Integral Equation 



A.l Constructing the Solution 

The Mellin transformation of a function f(p), see e.g. [46], is defined as 

oo oo 

M[f;s] := J dp/" 1 f(p) if J ^\f(p)\ 2 exists. 



(A.l) 



Applying this to both sides of (2.9) and using the faltung theorem [46, Chap. 4.8], one 
obtains the algebraic equation 



Qi.2 o \ 

M[t; s] = 8ttA M[KmK®(^- - 7 2 ; k,p); s] + 4=- 



-1)' M[Q l [x + -);s-l]M[t;s] 



(A.2) 

which is easily solved for A^[t;s]. Thus, one now only has to apply an inverse Mellin 
transformation, 



c+ioo 



*(p) = ^t / dsp- s M[t;s] 
Ziri J 



(A.3) 



where the inversion contour must be placed in the strip c in which all of the original Mellin 
transformations exist. 

However, there is no Mellin transform of /Cw in the limit 7, k <C p because it is propor- 
tional to (k) l /p l+2 . One therefore has to resort for the inhomogeneous term to a slightly 
more complicated, "half-plane" transformation [46, Chap. 8.5]: 



1 

Ik 2 f Ik 2 
AC[hm/C (z) (^- - 7 2 ;£;,p);s] := / dp p^ 1 lim/C (/) (-^ 7 2 ; k, 



P) oc 



k l 





00 



M + [\imK, { - l \^-- 1 2 ] k,p) ] s}:= /dp/- 1 lim /C W (^- - 7 2 ; k,p) oc 
fc-^o 4 J k— »o 4 



s-l-2 



s-l-2 



(A.4) 



These Mellin transforms A4_ and M+ exist for Re[s] > Re[l + 2] and Re[s] < Re[Z + 2], 
respectively. The solution to the integral equation is in this case given by [46, eq. (8.5.43)]: 



(T—+10O 



8ttA A4_[lim/C«(^ - 7 2 ;£;,p);s] 

ds p ; ; — ; 77 : + 



<T_ —IOO 

o^+ioo 



ds p 



cr_|_ —IOO 

ds p" 



8ttA A^ + [lim/C«(^ -7 2 ;fc,p);s] 

k^O 



S(p) 
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(A. 



where ct_ > Re[/ + 2] and a + < Re[Z + 2]. The denominator is simply the Mellin transform 
of the resolvent of the Faddeev equation. Not surprisingly, it determines the asymptotics of 
the solution. The function S(p), determined by the boundary-conditions, is in general an 
analytic function in the strip <r_ < Re[p] < a + in which the integration contour lies. 

We thus see that the particular solution is finally defined everywhere except at those 
points Re[s] = Re[Z] ± 2 where A4±[lim KM'] does not exist. It is not necessary to perform 

the contour-integrations leading to an analytic solution here. Rather, we note that 

'i(P) = £*-«W7 withl = -^- (-1)' M[Qi (x + - ) ; sf (A)] (A.6) 

~( p S l W+ 1 V37T V 

with some fixed coefficients q with which the zth zero Sj (A) of the denominator in (A. 5) 
enters at fixed (Z; A) - unfortunately, there is no closed form for these residues. We used that 
because Qi(x + 1/x) is real and symmetric under x — > 1/x, the zeroes in the denominator 
of (A. 5) come in quadruplets {±Sj (A); ±Sj (A)}. Only the := s closest to —1 is 
important for the amplitude at large p, as it provides the strongest UV-dependence. Notice 
again that only those solutions exist which do not diverge as p — > oo and for which the 
Mellin transformation M.[Qi [x + ~) ; s] exists as well. 



A. 2 How To Do an Integral 

To obtain the zeroes of the denominator - or equivalently equation (2.18) in the main 
text - we now perform the Mellin transformation of Qi [x + ~). First, one represents the 
Legendre polynomial by a hyperbolic function [27, eq. (8.820.2)], and then uses in turn the 
series-representation for 2 F\ [27, eq. (9.100)]: 

n ( l \ v^r[Z + l] / iv w r ^ + 2 Z + l . 3 / iV 2 i / a ^ 
Q \ X+ x) = 2mr[Z + |] { X+ x) ^ [ —>—> l+ 2-\ X+ x) ] (A - 7) 

v^r[z + i] y. r[| + i + n]r[^ + n] / lyC^+'+D 



2 ,+1 r[f + i]r[f] ^ r[z + | + n]r[n + i] 



This series is convergent because (x+ -) 2 < 1 for all x G [0; oo], cf. [27, eq. (9.102)]. Now, 
perform the Mellin transformation of each term using [27, eq. (3.251.2)]: 



This integral exists for Re[2n + Z + 1] > |Re[s]|, and hence for sufficiently large n, i.e. for 
an infinite, absolutely converging sequence. By analytic continuation, the result can thus 
be shown to be correct for all n. After a few simple manipulations also with the aid of the 
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doubling formula [27, eq. (8.335.1)], one can re-sum the series again: 



/- 9 -(Z+2) 1 I 2 2 J P 

2 fpUl 2 ^ 



/ + s + l Z-s + 1 2/ + 3 1 



Inserting this into (A. 6) leads to the algebraic equation (2.18) for the coefficients s/(A) given 
in the main text. 

When does a solution to the homogeneous version of (2.9) exist? In general, no arbitrary 
homogeneous terms can be added due to Fredholm's alternative: A non-zero solution exists 
for a given boundary condition either for the inhomogeneous or for the homogeneous integral 
equation. This follows also from the considerations leading to (A. 5) because the two regions 
in which A4-[hm /C^(^ ■y 2 ; k,p); s] and A4 + [lim /C w (^ 7 2 ; k,p); s] are defined do in 

k— >0 k— +0 

general not overlap, so that S(p) has no support. 

However, when M.[Qi (x + ^) ; s] itself exists, then the homogeneous version of the inte- 
gral equation has a solution. In that case, (A. 8) exists for each n, and in particular for n — 0, 
so that one must have Re[Z + 1] > |Re[s]|. As shown in Sect 2.3, the kernel is then singular, 
circumventing Fredholm's alternative. Danilov [10] discussed the case (Z — 0; A = 1), where 
the Mellin transformation is listed in [27, eq. (4.296.3)] with the constraint 1 > |Re[s; =0 ]|, 
consistent with our result. 

To summarise, all Mellin transformations of the particular solution in (A. 5) are well-defined 
for all (s, I, A) except for the driving term, and for the back-transformation (A. 5). This 
constrains the values of s to: 

Re[s) ^ Re[Z] ± 2 , Re[s) > -1 (A.10) 

The homogeneous part of the Faddeev equation has in general a solution only if the kernel 
is not compact. This is found for 

|Re[s]| < Re[Z+ 1] . (A.ll) 
A.3 The Full Off-shell Amplitude 

One obtains the solution to the full off-shell Faddeev equation (2.14) easily as follows: 
Replace 87rA M±[lim /C (/) (^ - 7 2 ; k,p); s] in App. A.l by 

k — ^0 

M[8nX "^T Ql (f + 1) ; s] = 87rA F ~ 2 M[Ql ( x + 1) ; s ~ 1] ' (A ' 12) 

The asymptotic form given in (2.15) follows now from the analogue to (A.5/A.6), keeping 
in mind that the contours can only be closed in the positive half-plane when k < p, and in 
the negative one when k > p. Notice that both off-shell momenta must obey the integral 
equations (2.14). 
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